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Abstract 


MontliXy mean simulations of the global atmosphere 
for October 1976 through Pebruary 1977 j initialized v/ith 
data from the beginning of each month, were computed with 
a new coarse-resolution general circulation model deTeloped 
at the Goddard Institute for Space Studies » Measured in 
terms of the monthly mean fields of sea-level pressure, 

850 mb temperature, and 500 mb height, the simulation 
skill of the model, which is still under development, was 
found to be inadeg^uate, thus far, compared with clima- 
tology. 

Substitution of observed monthly mean sea-surface 

temperatures (SSTs) as lower boundary conditions, in 

place of climatological SSTs, failed to improve the model 
* « 
simulations, I'llhile the impact of SSI anomalies on the 

model output is greater at sea level (v/here anomalously 

cold and warm water generate higher and lower sea-level 

pressures, respectively) than at upper levels (where it 

is negligible), the impact on the monthly mean simulations 

is not beneficial at any level » 

Shifts of one. and two days in initialization time 

produced small, but non- trivial, changes in the model-, 

generated monthly mean synoptic fields, Mo improvements 

in the mean simulations resulted from the use of either 

« 

time-averaged initial data or re-initialisation with 

\ 
\ 
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ijime-averagefi early model output. 

Qlie noise level of th.e model, as determined from a 
miiltiple initial state perturbation experiment, was found 
to be generally low, but with a noisier response to initial 
state errors in high latitudes than in the tropics. How- 
ever, the influence of random initial state errors on the 
monthly mean simulations is negligible compared with the 
large-scale simulation errors, indicating both stable 
model behavior and the need for further model improve- 
ment. 


% 
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Introduction 


A new global general circulation model (GCM) V7itli 
coarse lioriaontal resolution bas recently been developed 
at tbe Goddard Institute for Space Studies (GISS) (Hansen, 
1978). Designed primarily for climate investigations, 
tbe model is currently capable of generating one day of 
simulation on an 8 x 10 degree latitude -longitude grid in 
four minutes (on an IBM 360/95 computer), which is an 
order of magnitude’ faster than the ^ x 5 degree GISS 
model (Somerville et al., 1974) from which it was derived. 

Diagnostic tests of the new ''Climate Model" have been 
encouraging enough to warrant a preliminary evaluation of 
its capability as a prediction system, Oiie high speed 
of the model maJces it particularly attractive for long- 
range forecasting studies. Bor example, a 30-day fore- 
cast can now be executed in two hours. We have, there- 
fore, undertalcen to use the model for some monthly mean 
prediction experiments similar to those carried out vfith 
the GISS model. (Spar et al., 1976; Spar, 1977 a, b; Spar 
et al,, 1978; Spar and Luts, 1978), Because these fore- 
cast experiments are not conducted in an operational con- 
text, they are referred to here as simulations. 

Dike the GISS model, the climate model is a global, 
spherical coordinate, primitive equation system, divided 
vertically into nine dynamically active layers, with top 
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a,t XO iiib (p^)f Pottow at the earth’s aaT?face (surface 
pressure, p^), ajid slgm (^), the vertical pressiire cq-« 
ordinate, defined as (p - p^)/(pg ~ (iPor radiation, 

hut not djoaomical, calculations, three additional levels 
have been added in the stratosphere above 10 mh») Hori- 
aontally, the coarse-mesh version of the model used in 
these experiments employs a 24 x 56 gridpoint array sym- 
metrioal about the Equator, corresponding to intervals of 
approximately S*^ of latitude hy 10^ of longitude, mid a 
12-minute time step# 3?he Ara»ca\\ra (1972) computational 
scheme, including ’*5}ASU-Matsimo" extrapolation, was used ’ 
in these calculations, as in the GISS model (Somerville 
et al,, 1974)# (Eor this experiment, the computation 
time was eight minutes per simulated day* In more recent 
versions of the model, the running time has been cut in 
half with the introduction of leapfrog extrapolation and 
a 15-minute time step*) Because of the large grid sise, 
it is necessary to assign to each gridpoint fractional 
values of ocean, land, snovj, and ice representative of 
the areas surrounding tlie point for use in the physical 
calculations# 

Snow cover in the model is computed from the surface 
air and groimd temperature calculations, and surface 
.albedo is variable. Both solar and long-wave radiation 
calculations malce use of a generalisation of the 



^-distribution metixod" (Ijacis and Hansen, 1974), which 
includes multiple scattering and talces into account HgO, 
OOg, Og, trace gases, and aerosols* Glouds are treated 
as nQU’-blaok bodies with albedo dependent on zenith angle* 
Pry and moist convection in the model ai*e computed by a 
method based on the spatial variance oi* static energy and 
a form of coxiveotive adjustment, ^Dhe model computes pre- 
cipitation and ground wetness, and sub-surface as well as 
surface temperatures over land and ice. Sea-surface tem- 
peratures aj?e prescribed as boundary conditions* 

She purpose of the present study was to determine 
how accurately the coarse-mesh climate model simulates 
observed monthly mean states of the atmosphere from given 
initial and surface boundary conditions* It should be 
noted, however, that the model has been undergoing con*- 
tinuous development since these experiments v/ere started, 
•and that the results reported here are not necevssarily 
representative of its ultimate performance* 

Ihe period selected for this test was the anomalous 
winter of 1976-1977, which was characterised by unusually 
cold weather in tlie eastern United States with abnormally 
high temperatures in the v/est, I'our groups of eacperi- 
ments were carried out, all based on global data for the 
period October 1976 through February 1977 px*ovided by 
the national Center for Atmospheric Research (HOAR) and 
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tli.0 Goiiiiex' (BiO), and dej^lVGd 

,t3^om o^omtiomX MG anaXyaes^ In tlxe G:>spc53^imont» 

Xniijidl daiJn 00 GM' on tlio do^r oaoh moiiijh 

vfexo naod to XnltiaXlso ^1 '?q montli^Xong" fo5?ooast inma* 

At tHiQ ond of oaoh nun, tlio lS*-honnXy outputs of tlio model 
vono atens^^od to pnoduoo a p^^o^lotod montIil;ii^ moan state ^ 

In tlUs osjponimont, olimatolo^slQal moivldily avonago soa- 
sunfaoo tomponatxire (SB9?) ca\d soa^loo fields from tlio 
teid atlas (AloKmidon ai^d Mobley, 1974) were used as sur- 


face bcnmdany conditions. In a second e:;:periment, to 
moasiu'o tlio impact on tlio simulations of anomalous ocean 
temperatures, tbc five forecast rmis vjoro repeated with 
observed montlily moan SSO} values derived from satellite 
radiometer moasuroments (l^rowor ct al,, 197^) used in 
plaoo of tlic olimatologiocOI. SSI's* 

A tlxird group of experiments was carried out to 

1 

determine the offoet on tlie montlily mean simulations of 
a systematio alteration in tlio iiiitial synoptic pattern. 


such as n\ight result from a shift in initialisation time 
or *Uio smootliing of initial conditions, IVo of tlie montlily 


forecast runs (for October 1976 and danuaey 1977) wore re- 


peated twice, starting with observed initial data for 
00 GMf on the sooond and tliird day of oaoh month, respec- 
tively, rather 'bJian ’Wie first, Avera^;es for correspond- 
ing periods to the end of tiio month were tlion compared in 
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order to assess the sensitiTity of the model-generated 
monthly means to the arhitrary choice of initialization 
time* In a related experiment^ the data were areraged 
over the first 2^ days of the months to oDtain a set of 
time-averaged initial conditions* !Hb.e forecast , computa- 
tions v;ere then repeated with these as .initial fields to 
determine whether such smoothing of the initial data -vfould 
result *n more realistic monthly mean simulations than 
those generated from the “noisier" instantaneous initial- 
izations* Another similar experiment was also conducted, 
in vAiich the model was re-initialized with time-averaged 
output from the first five days of the original forecast 
runs before the completion of the monthly simulations* 

A fourth set of computations was performed to mea- 
sure the "noise level" of the model (Ohervin and Schneider 
1976 a, h)* In this experiment, which closely resembles 
the random initial state perturbation experiment carried 
out with the G-ISS model (Spar et al,, 1978), Initial con- 
ditions for October 1976 were contaminated foin? times 
with different random error distributions, resulting in 
a total of five simulations for that month, Ih.e dis- 
persion of the simulations represents the inherent un- 
certainty, or "noise", of the model-generated monthly 
.mean fields associated with unavoidable random errors in 
the initial conditions. Such noise level statistics 
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5.ndicate not only tKe irreduci'ble miniranm error of the 
montlily mean atmvaations# but also the minimum deteotable 
signal from a prescribed climate change experiment with 
the model (CherYln aiicL Schneider, 19T6 b). 

Karee q.uantities were selected for analysis and 
evaluation; sea-level press^^re, S50 mb temperatiare, and 
500 mb height. Simulated and observed monthly meaaa fields 
of thC'Se variables were computed over the globe, but the 
main empliasis of the evaluation was placed on the Northern 
Hemisphere, with particular attention to the quality of 
the simulations over tlie North American quadrant during 
this anomalous winter, 31ie simulation skill of the model 
was evaluated numerically over seven regions, ranging in 
siae from the United States to the whole globe, in terms 
of area-weighted root-mean-square (rms) errors and SI » 
(gradient) skill scores (QJeweles and Wobus, 1954), 

To provide a standard for evaluation of the model 
simulations, monthly climatologioal fields of the three 
variables, interpolated from data fmmished by NCAR^, 
were also evaluated as ’’forecasts” of the five observed 
monthly mean states. Comparison of the error statistics 
of the model simulations with those for climatology pro- 
vides an objective measure of the model's capability of 

reproduoing ollmatic anomalies, 

5 

Ihe NGAR climatology is derived from a variety of 
sources as described in Crutcher and Meserve (1970) 
and Jeime et al, (1974- )* 


8 



I?asiG SimuXation 3Sx :’nGriiften1;^ ‘ 

ymiiii II 1WI BiiwT imiiiiiiTMiim ■ T II II ■ II i lauM minn-B 

S}liG iJii’feiaii.EiGd wi til gloTial data 

interpolated :£;eom tliG op erationai MHO analyses for 
00 0143} on tile first day of eaolx of tixe five months ^ 
Ootohor 1976 tlirougli S'e’bruary 1277* Sea-snrfaoe tem- 
peratures and sea-ioe looationa, hased on tlie monthly 
climatological data of Alexander and Hohley (1974) » hut 
interpolated daily, were specified as surface, boxmdary 
conditions o 

l‘he chnractei'istics of the anomalous cold vdnter of 
1276-1977 began to appear over North Amearica during the 
latter half of September (iPaubensee, 1S7S) in the form 
of a deep tropospherio trough near the east coast of 
North America with a strong ridge in the west. As a i'e- 
suit of this upper air flow pattern, which persisted 
tlirough January, record lov/ temperatures developed .over 
the eastern United States, vdiile tempei'atures were ab- 
normally worm in the west (Wagner, 1977 a; Dickson, 

1977 aj SCaubensec, 1977 j Vfegner, 1977 b)* In Januar^r 
the wave pattern amplified still fiirther (Wagner, 1977 b), 
leading to repeated advection of Arctic air southward 
into the eastern United States, wlxere surface tempera- 
tures fell far below normal t notably in the Ohio Valiev^ 

^Wapted in part from a master »s thesis sxibmitted to the 
Gity College by Hobert d* liutc« 
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I^CQduci^ig one of the coldest months on necond* 

- I'ehnuEiny h^jonght. aho modification and eyentnally’ 
a hrehicdom of this anomalous situation (Dickson, 1977 h)* 
3}he mean nidge oYen the Pacific ITQnth’,^est moved east\?and 
and the mYe pattenn flattened* Pantioulaxly during the 
last week in Pehruany, the mid-tropo spheric circulation 
OYer ITorth ihaerica changed abruptly. She large-amplitude 
VaYe pattern, which had dominated the winter season, 
hreke dom and gave way to fast westerlies with fast 
moving storm systems traveling across the coimtry. ¥lth 
westerly flow now dominant, \'^arm air spread quickly across 
the country and temperatures rose above normal in the 
eastern United States . 

How well the model simulation reproduced the anomalous 
circulaticn and temperature fields diiring the winter of 
1976-1977 may be seen in Pigures 1, 2, and 5, which show 
the distribution over tlxe northwest quadrant of the earth 
of sea-level pressures, S50 mb temperatures, and 500 mb 
geopotential heights plotted in units of mb - 1000, 
degrees Celsius, and deoameters - 500, respectively. Por 
each month the observed field is shomi at the top and the 
simiilation at the bottom of the figure* 

In Figures 1, 2, and 5, North America lies in the 
.center of the region, with the North Atlantic Ocean in 
the right-hand third of the map and the eastern North 


4 


■ a 

■ : 3-. ■■ 

'■4,- 


" * 

' ," T" 

«■ 

■:«■ ■■ ■ 

10 

10 

B» 

&» 

I® 

'S3' 

D« 


a® 

a ® 

8< 

84 

13 

i,e 

. a* ,' 

t7 


Mr* 

I® 

■ 

— ,14 — 

•34 " 

-tor- 

33 

ifrj 

•<0 '• 

s® 

✓ 

1® 

!*■ 

83 

1®' 

8* 

84 

84 


OEIGIHAL PAG® IS 
0® POOB QUALITY 


f sv 

K . 


^ w 4, 

II „!?- 


» } • 4 . _ - — ~ - 

ai II ‘'*.•11 * » ss 8* la* 83 s« *• 1 * * * i a ^ T ••s ■*-*.„* j 

-- ■ \ I ' 


«* «• 4 

H 4 S * , '1^ 8* 1* *T ^io * ^ 

je . i 4 '**‘Ta 1 * Vsy'^i# /fS\ la iw w i® tr \ i» *» jJif‘'"'8» .'n 

. ■-«" ^ ' j ^ \ ” *** 

f» * 0 ^ i® t® *7 I* fir s« 19 a® Jtt/*' IT 1 ® i' 33 taV. '“y. 

*“ U '* V ^ V V )/ 

IS i« 1® t» 1® la *« ^-44 V*^— IB- .tui^ *3 IT 19 \*|- *x £©//l5 

Jg ts I« 13 I* as is tsj I* la- 13 is as'’v it w l* 

_ ->1 - »* *•* - _ 


-V ^ 

£ I 9 /»4 IS IS as 13 / it 10 

^ V- «— ♦ ■"** ^ \ ^ . 

■ .-^E* 'If 89 1* IT tOv. 1*^ —IR. 


V. ^ I S’" 

•x ^ ^ -3 ■>. „-Z. 

■V .. — — 


J S a 3 4 • 0 T • * 10 It 13 

at »2 a* ea ** *a aa £z «z ax aa *a s* 

cj al******? ”'"'£ij*'>^ as'^'^'-fctr- '"as t* as J ®i 

— 

*r fa i as ao 36 Jt9 st* SSPio' 


!2 *r ”***28 ^£*'***33 "i 32 JO 30 30 3* ' ** } ** 

ft as ,-'»''«-~rB''N3[8^-''''3Ey''«s st 3^ at*^^ — \ ] 0 ^ *o 


/ioN 


u » «- ae to to to is so 

IS 9 « »« II 11 


is T 
••too 


9 « to II II IS It 

•140 •‘**0 -10® -80 


to 

\ T 

a 


“ 

It 


IS 

2> 

1® 

14 

AX 

y* 

14 ^•'12 

* 4^ 
\ 

\ 

£1 

a« 


1® 

11 

.-to- 

SI 


iAi-» Me«v 

(. 

«SaO li. 

mi*t s» 

T«.s i : 
04 uS ' 

«8«T f 

SQaO 1 

43<3 1; 

3S«X !■: 

»r,n- tr 

19«* 1* 


•to 

11 


13 

13 

14 

13 

IX 

tl.T 

i: 



S 








11 

10 

3« 

1^. 

13 

13 

13 

13 

^•O . 

U 









LAV 

»EA 

•*0 


-40 


»40 


-ZO 




11 

IS 

IS 

14 

t® 

a® 

IT 

t® 

lAf. 

H^Af. 

Z2 

S£ 

K£ 

*S 

gz 

xz 

XX 

XX 

. 90.® 

2:^ 

<£j 

St 

•-a 

/«A 

• i 

SZ 

ztt 

y “X 
yZ7 

»-2J/ 

l/' f~ 

t.zv 

flZ«£ 

SA 


Tt<3 27 

AS«S £0 

se.r 13 

•EOaO as 

*3.3 is 

33.Z 14 

aTtt* S3 

IOa« 11 


Figure !• 


Mouthiy mean ’-sea’-' level pressures ^ over the northwest quadrant. 
Top; observed? bottom; simulated. 

Units; mb - 1000. 4 mb isobars.' 
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Figure 1. (b) November 1976. 
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Figure 1. {c) December, 1976, 
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Figiire 2. (b) November 1976. 
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Figure 2. (c) December 1976 
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Figure 2, (e) February 1977. 
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(c) December 197-6 


A .. 

t 

: : A ■ 


*6 , ■ ■ 

4 

. .4 , 

- V : 

m.' 


9® 

US .. 

la 

13 

14 

*» 

as . 

*y 

la 

14 

as 


' 3l'_ 


•It 

Si 

*L, 


«*- 

31 

4“L, 



■3 » : 

wiU, 

3| 

at 

XI 

>3 

a« V 

34 

33 

Jy-"' 

a« 

a* ' 

Mir 

34 

xa . 

B4 

24 ' 

aa 

. aa ■ .. 

aa 

'£3 

ax 

- 1 
« t 

J 

«J 


«* V n 

at 4a 

ao 

tt 

is 


£r a* . 4« ax 


x» 


r' 


tS'< *a V ay 

) «/ 

# / y / 


■U 33 
A4 


4T 

N. 


I* 


1* 


a® i. aa \,i* ise ja 

^fik. ' t 


•'s. 


,-S. 


ta 


4® 'via 

/ - ^ ''v^ 'v 

sa »3 V 40 'S.** '^.»T 


t» 

»r 

ta 

ta 

^ V* 

a« 


I*'! at 

y 


XB aa 


!xl <3 

ar 

as 

1* Nat 

*• j 

f 33 

V 


( 

1* ix' 

^-.4w«r^'^- ^ 

jy 







aa aa 
»? a* 
34 aa 
•Xi' <4 
»T or 


aa 

as 

“aa 

42 

34 



la a a ajg aw _ 


24 ^^ 0 — ST aa **‘“”-'4?^ ** 

''•S ^ 

'3fi— ar aA^^—35 


s,w* Kt 

'Mtj® 

aa«a 

V««a 

4A.S 

8B.T 

so.® 

43*3 


1* 


4* 

,-*•***"* 

_^«7 

****^ 

os'- 

•kO^Wav. 

-*< 

1— ^18-.*— 



*3 

**'**.. 

3S.X 

14; 


j4>. ■ 

7* 

>■ •,* 

74 


74' 

7* 

73 

• 4^ 1 

73 

TT 71^ 

«3L4 •-* -•- 

.1? 

^ .-*4- — 


77 7« 

£7^4 

la 

«3 

es 43 

«X 




42 

«« 

•4 

04 

S7 

BT 

*3 

84 

«3 

03 S3 

i®.4 

14 

se 

»s sr 

•4 

84 


•t 

89 

«4 

87 

•r 

8T 

87 

84 

'84 

84 

84 ar 

11.7 

la 

a® 

•7 87 

•« 

84 

«« 

as 

84 

84 

84 

<& 

«34 

84 

84 

84 

04 

87 87 

3.® 


«ieo 

-140 


*140 


-izo 


too 


*80 


•40 


*40 


*20 

tAT 

J 

t 

a 3 

«l 

a 

4 

y 

e 

9 

IP 

11 

* 12 

13 

14 

18 

18 

»T 1« 

LAT. 


*0i>4 

OXafi 

r*.a 

4«aS 

a«,r 

so*® 

43.3 












*— s*«« 

» 

*— •■■ 



, 


IT 

^4®. 


37 

83 

34 

" S2 

~"m£' 


<6 

43 

4A 

y.^-a3r"" 

or 

^5H_ 

__S* 

DS 

34 

3S.X 














* _ 







14 

44 «3 

43 43 

43 

43"^ 

“40“ 

40' — 

-fiff* ** "e* '*43 

4S 

“ ** ■ ■■»■«■ 4 

aa 

44 

44 

44 

83 

43 

27.4 

19 

■ 73 

tjT 

'**'70—' 

Vl 


*—-7,t — '“jr" 

"* 72 

74 

7S 

73 

7* 

72 

72 

~7t~' 


’ "rs" 

73 

I9«4 

H 

78 

TT 

74 

74 

74 

77 

78 

TT 

78 

78 

77 

77 

74 

74 

73 

74 

77 

T* 

ll.T 

13 

7* 

. T4 

7* 

78 

78 

78 

78 

78 

78 

rr 

77 

77 

78 

77 

77 

78 

78 

78 

3.9 

. 

•>t«0 


-140 


*140 


-120 

IN 

■108 


“00 


•"40 


-40 


*20 


lAT 


Figure 3* (d) January 1977. 
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Figure 3. (e) February 1977 


Pacific Ocean on the left, longitudes are labeled at 
tlie bottom of each map, V7hich extends from 10° \i on the 
right to 180° ¥ on the left, the negative numbers denoting 
\iest longitude, fhe latitudes of the gridpoints are 
shown in the second column on the right of each aiap ixi 
degrees ¥orth. (Kote that the latitxide interval is ap- 
proximately, but not exactly, eight degrees.) Ihe first 
colniim on the right of each map shov/s the aonal mean 
value of the variable displayed within the quadrant, 

(G-rid numbers also appear, at the top and in the left 
column, ) 

In -Figure 1, sea-level isobars have been di’avm at 
an interval of 4 mb. It is apparent from an inspection 
of tlie maps that the model fails to simulate adequately 
the monthly mean sea-level pressure fields. Although ' 
the Icelandic and Aleutian lows, as well as the svib- 
tropical high pressure cells over the oceans and the 
North American continental high, are all reflected to 
some degree in the model smnulations, the quantitative 
agreement betx'xeen the observed and predicted mean fields 
must be characterised as poor. Most notable are the 
failures to simulate the abnormally deep iUeutian lows 
in January and February, and the strong pressu-re gradients 
in the North Atlantic in October, November, January axxd 
Febi’uary. In general, the model-generated monthly moan 
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eea-leTel pressures are too high in high latitudes and 
too low in low latitudes* 

MiG 850 mb isotherms in Figure 2 are dravm for an 
interval of 5^^ 0. It is immediately apparent that the 
predicted temperatures are generally too low, especiaH.ly 
in high latitudes. The cold Arctic may be due to the 
inadequacy of the eddy transport of sensible heat, a 
problem that is also present in the G-ISS model (Stone, 
et al.j 1975). However, the present model simulation 
is too cold in the tropics as well, which suggests some 
other defect, 

The dominant characteristic of the observed tempera- 
ture field over the northv;est quadrant from October 
through February was the large contrast betv/een high 
temperatures over the west coast of Horth America (e.g. , 
120^ W) and low temperatures in the east (e.g., 80° V/), 
which became most pronounced in December and January and 
began to amelioi'ate in .February. The model docs simulate 
some aspects of the cold wave in the east, notably in 
November. However, the December and January simulations 
are less satisfactory, and in February, when the cold in 
the east had abated somewhat, the model predicted the most 
severe conditions. Vftiile the model does depict large 
negative deviations from normal, it fails to simulate the 
phase opposition between the cold east and warm west, and 


in Eetouary it grossly exaggerates t’ne cold anomaly in 
tile east. 

aHie observed and predicted 500 mb iieigbts, together 
with msimally-dravni 100 m contonrs, are displayed in 
figure 3- iUhe observed fields shovf the persistent sta- 
tionary wave pattern (western ridge, eastern trough) men- 
tioned earlier. She positions of the ridge and trough 
over IJorth America are reasonably well-simulated in the 
prognostic maps. However, the amplitudes of the mean 
monthly waves, especially in Hovember, December, and 
January, are not adequately reproduced, and the predicted 
flow in those montlis is much more zonal than the observed, 
nihe February 500 mb simulation is much closer to the ob- 
served field over Horth America. In all five months the 
predicted 500 mb heights are generally too lov; compared 
v;ith the observed values, which is consistent with the 
fact that the troposphere is too cold in the model simu- 
lation. 

For a q_uantitative evaluation of the model output, 
simulation error statistics v;ere computed for seven 
regions of the earth, ranging in size from the United 
States to the total globe and defined as follows: 

(1) United States (27° H to 51° H and 130° W to 70° V/) : 

(2) East Pacific and United States (27° N to 51° H and 

180° U to 70° (3) Horth America (27° M to 74° K and 
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130° W to 70° W); (4) Europe (35° K to 74° N and 10° W 
to 40° E)s (5) iPropios (20° to 20° S bver all longi» 
tudes); (6) Northern Hemispliere (4° H to 82° N over all 
longitudes) j (?) G-lobe (82° S to 82° E over all longitudes). 

{Qlie bias of the model is indicated in liable 1, v/bich 
lists the area-weighted algebraic mean errors of the 
monthly mean simulations over the Northern Hemisphere, 

Also sho\mi for comparisons ai’e the corresponding values 
for a ’’forecast’' of climatology. As noted earlier, the 
simulated sea-level pressiu'*es in the Icelandic and 
Aleutian region were generally too high. Hov/ever, these 
local effects are overcompensated by underpredictioiis 
elsew’here, v;ith the result that, for the Northern lieni- 
sphei'e as a whole, the average bias of the soa-lcvcl 
pressure simulations is - 1.6 mb for the five months. » 

The 850 mb temperatures in the simulations are too lovj 
everyv/here, including the tropics, v/ith an average cold 
bias of - 3.5° C for the five months over the heiiii sphere. 
This result is hydrostatically consistent with the nega- 
tive bias of the simulated 500 mb heights, which average 
95 m less than the observed heights ovex* the Northern 
Hemisphere for the five-month period. Compared with 
that of the simulations, the bias of the KCAR climatology 
in Table 1 is negligible. This indicates not only that 
the hemispheric mean state v/as close to normal for the 
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gable .1* Algebraic mean difference (bias) between 

monthly mean fields over the Northern Hemi- 
sphere tot winter 1976-1977 • M: Model 

simulation - observed. G; .Climatology - 
observed. 


October November December January gebrtiary Average 

Sea-level pressure (mb) 


M 

- 

1.4 

- 

1.7 


1.3 - 2.2 

1.5 

- 1.6 

C 


0,2 


0.0 

- 

0.0 - 0.7 

0.3 

- 0.2 





CD 

o 

mb 

temperature {°C) 



M 


3*6 

- 

3.5 

- 

3.4 - 3.5 

3.7 

- 3.5 

C 

- 

0.4 

- 

0.7 

- 

0.4 - 0.5 

0.7 

- 0,5 





o 

o 

mb 

height (m) 



M 

- 

94 

- 

91 

M \ 

B6 - 100 

93 

- 95 

C 


4 

+ 

5 

+ 

4 +2 

5 

H- 2 
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period (see also, Atigell and Korshover, 1978), tut also 
tliai: the model exhibits a negative bias in pressures, 
temperatures, and geopotential heights relative to clima- 
tology as well as to the observed monthly mean fields. 

0?he rms errors and S1 skill scores for the five monthly 

k 

mean simulations, together wii:h the corresponding scores 
for the ’’climatology forecast”, are sho^TO for the seven 
regions defined above in (Tables 2 and 3. (The S1 score, 
a conventional measure of the difference between predicted 
and observed gradients, is a dimensionless quantity , v/ith 
a range from zero, for a perfect forecast, to a maximum 
of 200. Based on experience at IMG (Shuman and Hover- 
male, 1968), it is generally agreed that 31 scores less 
than 20 indicate virtually perfect synoptic agreement, 
v/hile scores greater than 70 represent "v;orthless” fore- 
casts, at least for 500 mb patterns.) 

Brora Tables 2 and 3 it is apparent that the model 
simulations are, in general, inferior to climatology. 

Only in the few regions and months indicated by asterisks 
is the error score of the simulation smaller than that 
of climatology. The rms errors, expecially of 850 mh 
temperatures and 500 mb heights, reflect, among other 
things, the negative bias of the model noted earlier. 

The large S1 simulation scores for sea-level pressure ^ 
indicate no skill, even where (as in Europe) they are 
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U) 


ggable Zo Root“meaU“sq_tiare (rms) eTjcoTs. of monthly mean 
model simulations (M)- and climatology (C) for 
winter 1976-1977. (See tesct for definitions 
of regions . ) 





Oct, 

Nov. 

Dec. 

Jan. 

Peh. 


Sea-level 

ores sure 

(gib) 




1. 

United States 


5.5 

5.1 

4.5 

4.5 

6,0 . 



C 

1.6 

5.1 

2.2 

4-0 

5.1 

2. 

E. Pacific - U. S. 

n 

4.5 

7.2 

6.8 

10.5 

8.7 



C 

1.9 

5.8 

5.4 

7.1 

6.1 

P • 

north America 

M 

6.4 

7.1 

7.4 

8.9 

7.5 



C 

2.0 

5.1 

2.1 

5.8 

5.5 

4. 

Europe 

M 

5.5^ 

Q.l 

4.4* 

5.1* 

7.5 



C 

6.8 

2.0 

4.7 

5.4 

4.9 

5. 

Hiropics 

M 

5.5 

5.7 

2.8 

2,9 

4.1 



0 

1.5 

1.7 

1.8 

1.2 

1.6 

6. 

Northern Hemisphere 

M 

5.5 

7.2 

6.1 

7.2 

7.8 



C 

5.0 

2.9 

5.5 

5.9 

4.4 

7. 

Globe 

K 

8.5 

6,7 

7.0 

7.4 

8.2 



C 

5.0 

5-2 

4.6 

4.7 

5.7 


* denotes simulation values 

smaller 

than 

those 



for climatology. 
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gable 2 , (Pont ' d, ) 





Oct. 

Hoy, 

Dec. 

Jan. 

Eeb. 


850 mb temperattire 

(°Q) 




1. 

United States 

M 

3*7 

4.4 

3.8 

5.3 

6.1 



G 

2.7 

2.8 

2.4 

3.2 

2.4 

2; 

E. Pacific - U, S* 

H 

3.7 

4.0 

3.4 

3.8 

5.5 


‘ 

G 

5.2 

2.9 

2.3 

2.8 

2,2 

^ • 

Uorth America 

M 

6.8 

6.7 

5.9 

6.2 

7.1 



C 

2.4 

3-0 

2,5 

3.9 

3.2 

4. 

Europe 

M 

5.5 

5.1 

4.0 

5.7 

6.6 



0 

1.9 

1.4 

1.7 

1.7 

2,7 

5. 

Tropics 

M 

3.6 

3.2 

2.9 

3.4 

5.6 



C 

1.8 

1.6 

1.2 

1.4 

1.7 

6 * 

Northern Hemisphere 

M 

5.1 

5.2 

5.5 

5.9 

5.4 



0 

2.2 

2.1 

1.9 

2.7 

ro 

• 

7, 

Globe 

M 

6.8 

6.3 

6.2 

6.5 

6.6 



C 

2.7 

2.1 

2.0 

2.4 

2.2 
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gable 2\> (Oont^d., ) 





Octo 

Ho-v, 

Dec. 

Ian. 

Eeb, 


500 mb height (m) 





1. 

United States 

W 

109 

115 

105 

133 

149 



c 

48 

74 

57 

85 

49 

2. 

E. Pacific - U. So 

M 

106 

108 

102 

154 

136 



0 

46 

66 

54 

101 

74 


ETorth imerica 

H 

119 

124 

109 

146 

139 



0 

44 

71 

58 

89 

51 

4. 

Europe 

M 

157 

160 

115 

160 

162 



0 

62 

. 27 

60 

58 

63 

5. 

Tropics 

M 

85 

82 

78 

86 

90 



0 

, 16 

16 

15 

17 

20 

6. 

ITortliern Hemisphere 

M 

104 

114 

112 

131 

. 119 



C 

59 

40 

43 

75 

55 

7- 

(rlobe 

11 

102 

107 

99 

114 ■ 

109 



0 

49 

47 

56 

59 

44 
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gable SI skill scoi'es ior monislaly mean moae'l sirnu- 
lations (M) and climatology (C) for v/inter 
197S-1977. (See text for definitions of 
regions • ) 

pot?^ 3j?0Vo Sec* Jan* Feb* 
Sea*-level •pressure 


1, United States Ti 

0 

2. 1* Pacific - U* S. M 

0 

5. North. America M 

• U 

4* Nurope M 

0 

5 , Iropics H 

0 

6, Northern Hemisphere M 

G 
M 
C 


95 

120 

115 

156 

105 

47 

77 

52 ‘ 

90 

94 

112 

140 

155 

156 

, 159 

51 

78 

51 

SI 

97 

106 

.114 

115 

158 

91 

58 

82 

54 

95 

95 

85* 

125 

85* 

90 

108* 

107 

40 

101 

78 

111 

72 

69 

72 

68 

70 

54 

54 

57 

44 ' 

49 

95 

102 

90 

89 

87 

65 

55 

51 

64 

62 

88 

85 

39 

84 

85 

55 

52 

65 

56 

52 


7 


G*lobe 



gable 3. (aoirb^d.) 

^wn ffr "■ ~i • • 'iflirinr i ’ ni nr., ' tt l^■^il^ 




Dot. 

Hqv. 

3)0Q* 

Jan* 

Peb 

S50 mb 

teitmei^aiiui^G 





TJnited States 

M 

48 

37^ 

38 

36* 

44 


0 

A3 

46 

32 

B5 

34 

1* 3aoi^lc « b* S* 

M 

41^ 

50^ 

34 

42* 

Al 


0 

43 

42 

30 

49 

36 

bo3?tli Amorioa 

M 

S3 

30 

38 ' 

53^'^ 

61 


0 

41 

40 

40 

37 

4-0 

Biu^ope 

K 

61 

6? 

76 

95 

65 


0 

46 

* 34 

46 

49 

59 

gx'opics 

H 

73 

S6 

S3 

9S' 

97 


C 

5B 

64 

62 

70 

67 

Kortliern HaniispberG M 

66 

63 

67 

6? 

62 


0 

49 

47 

59 

31 

43 

Crlobcs 

M 

61 

61 

62 

64 

61 


a 

40 

41 

37 

44 

58 


ai 


galjle (Oont^d.) 





dot. 

Bov, 

Dec. 

Jan. 

Beb. 


500 mb liei.^ht 






1. 

United States 

. M . 

39* 

45* 

38 

50* 

50* 



c 

41 

54 

32 

52 

36 

2. 

B. Pacific - U. S- 

M 

45 

51 

43 

65 

43 



0 

37 

45 

51 

53 

40 


JTortli Jimerica 

-M , 

41 

49* 

49 ' 

65 

42 



C 

57 

53 

41 

58 

57 

■4* • 

Europe 

M 

65 

76 

54* 

71 

46* 



C 

57 

54 

66 

67 

61 

5. 

ICropics 

n 

55-x- 

71* 

S3 

72* 

102 



0 

66 

75 

72 

79 

S3 

6. 

Uortliern Hemispliere 

M 

52 

71 

67 

67 

. 62 

« 



G 

47 

42 

43 

57 

48 

7. 

Globe 

M 

45 

50 

52 

54 * 

49 



C 

38 

37 

43 

45 

58 
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smallea? tixose for climatoXogy* Over i;lie United Si;ates, 
t!ie S1 simnlation scores for tile 850 inU temperature and 
500 mb lieiglit fields arsi on tl^e ayerage, sliglitly araaller 
tlian tliose for climatology (40 vs, 4& and 40 vs, 45? re« 
speotively), indicating some skill in tlie reproduction of 
those patterns, possiUly due to persistence of initial 
conditions during this aiion^alous winter. However, there 
is no evidence of this skill over the northern Hemisphere, ■ 
for which the average S1 simulation score for 500 nib 
heights is an unsatisfactory 64 compared with al^ average 
of 47 for climatology. 

In general, it must be concluded that, despite the 
model^s realistic simulation of the large-scale meridional 
structure of the atmosphere, the monthly mean simulations 
do not adeg.\iately reproduce the observed synoptic fields, 
?Ehey are, in fact, inferior to climatology, ind3. eating 
that, at this stage, the coarse-’mesh model is not yet 
capable of duplicatiJig realistically significal^t depai’- 
tures from climatology of monthly mean synoptic patterns, 
‘JHais conclusion may be tempered by ,the fact that clima- 
tological SSfs rather than observed sea temperatixres were 
used for the calc\\lations of the -vertical fliuces of heat'' 
and water vapor over the ocean sixrfaoe in the simulations 
, above , 


H5 


Impact of Sea-surf aoe. temperature (SSg) Anomalies 

Real data experiments with the b-ISS model (Spar and 
Atlas, 1975; Spar et al., 1976; Spar’ and huts, 1978) have 
thus far failed to reveal any K.onsistent heneficial im- 
pact of observed SSt anomalies on model atmospheric simu- 
lations* Nevertheless, the viev is widely held (see, 
e.g., Harnaclc and Landsberg, 197S) that air-sea inter- 
actions play a signif icant role in the generation of at- 
mospheric anomalies on monthly, seasonal, and even annual 
time scales, and that persistent SSt anomalies may force 
the atmosphere into anomalous circuJLation and temperature 
patterns. In a recent review of the North American ab- 
normal winter of 1976-1977, Namias (1978) has again pre- 
sented interesting synoptic evidence in support of this 
concept. 

As a further test of the SSI hypothesis, the five 

simulation runs with the climate model were repeated vath 

observed monthly mean SSIs replacing the climatological 

values. The observed means for each month were computed 

from daily SSTs, which are derived from satellite radi- 

ometer measurements and provided by NOAA-^ on a latitude- 

longitude grid smaller than 1^ square (Brower et al., 

1976), Bor use in the model surface flux calculations, 

5 

'^National Oceanic and Atmospheric Administration, 
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4-poa.niJ averages of these fine-mesh monthly mean SSTs were 
determined for each coarse-mesh model gridpoint. 

The mean SST ‘anomaly field for each month is dis- 
played in digital form on the global maps in ]?igure 4 
(a-e), where negative (positive) numbers denote colder 
(warmer) than normal v/ater, in degrees Celsius. Blank 
areas on the maps represent the continents. Longitudes 
are labelled at the top and bottom, negative for west and 
positive for east, with the Western Hemisphere on the 
left. Latitudes are indicated in the second column from 
the right, positive for north and negative for south. 

The first column on the right lists the zonal mean anomaly 
for each grid latitude, (Also shovjn on the maps are the 
index numbers of the 24 x 56 gridpoint array.) 

The coarse-mesh climatological SST fields used for 
the calculation of the anomalies in Pigure 4 (and also in 
the basic simulation experiment) were computed by averag- 
ing the 1° grid data from Alexander and Mobley (1974) over 
the 8 X 10 degree box surroimding each gridpoint, and are, 
therefore, somewhat smoother than the observed SST fields. 
Also, the, observed SST tapes were' incomplete, with as few 
as 14 days of data available for the poorest month (Peb- 
ruary 1977). These factors may have contributed some 
irregularity to the SST anomaly patterns, although there 
is no reason to believe that the data presented in Pigure 
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Figure 4, Monthly mean sea-surface temperature anomalies {observed. — climatology) over the globe 
in degress C* 


(a) ■ October 1976 
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Figure 4. {b) Hovenbsr 1976. 
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Figure 4. (e) February 1977 


4 are incorrect, Hevertlieless, the SSI anomaly fields 
are relatively smooth, with certain persistent and co- 
herent large-scale patterns discernible in most months. 
Bins, cold water (negative anomalies greatex* than 2® C) 
is found almost every month in the tropical Pacific, in 
high latitudes of the ITorth Pacific, and in the G-nlf of 
CwUinea, while warm water (positive anomalies greater 
than 2*^ C) appears mainly in the western North Atlantic 
and in the Sea of Japan, Over most of the earth during 
this five-month period, the SSI anomalies were not 
greater than 2 ^ C. 

!nhe impact of SSI anomalies on the model-generated 
monthly mean fields ma 3 ’" be evaluated by comparing the 
two parallel simulations, based on climatological and ob- 
served SSIs respectively, for each month, Ihe c[Uantita- 
tive results of this comparison are shoim in lables 4 and 

5 for the Northern Hemisphere only. In lable 4 the rms 
errors and SI skill scores are compared for the two simu- 
lations, A, computed with the oJserved (anomalous) SSI 
field, and K (taken from lables 2 and 5), computed with 
the climatological SSIs, It is apparent from Table 4 that 
the use of observed monthly mean SSIs did not result in 
any improvement in the simulations over the Northern 
Hemisphere, and in many cases produced even poorer agree- 
ment with nature. Similar results (not sho\Mn) were found 
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feble A-, Rooij-meaw-sq.’aare (xms) errors ard S1 skill 
scores of monthly mean simulations computed 
with observed sea-surface temperatures (A) com- 
pared with errors (from lables 2 and 5) of 
simulations computed with climatological SSSDs 
(H). October 1976 tlirough February 1977. 


northern Hemisphere * 

1976 


1977 



Oct. 

ifov. 

bee. 

Jan. 

Teb7 

rms errors 

Sea-level pressirce (mb) 

A 

5.5 

7.4 

6.6 

7.4 

8.2 


Vi 

5o 

7.2 

6.1 

7.2 

7.8 

S50 mb temperature (°C) 

A 

5.2 

5.5 

5.5 

6,0 

5-6 


M 

5.1 

5.2 

5.5 

5.9 

5.4 

500 mb height (m) 

A 

105 

112 

112 

155 

116 



104 

114 

112 

151 

119 

S*l skill scores 

Sea-level pressure 

A 

98 

104 

92 

92 

95 


M 

95 

102 

90 

89 

87 

850 mb temperature 

A 

72 

69 

69 

75 

68 


M 

66 

65 

67 

67 

62 

500 mb height 

A 

55 

71 

70 

76 

61 

’ 

M 

52 

71 

67 

67 

62 
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laljle 5* Root-mean--sg.uare (rms) differences and S1 com- 
parison scores “between monthly mean simulations 
computed with observed (A) versus climatological 
(M) sea-surface temperatures. October 1976 
through February 1977 • Northern Hemisphere only* 





1976 


1977 



Oct. 

Nov. 

Dec. 

Jan. 

Feb. 

rms difference 







Sea-level pressure 

(mb) 

1.6 

1.7 

2,4 

2,0 

2,1 

850 mb temperature 

(°c) 

1.6 

1.5 

2.1 

1.8 

1.7 

500 mb height (m) 


12 

17 

24 

24 

21 

S1 comnarison score 



• 


• 


Sea-level presstire 


55 

58 

58 

41 

55 

850 mb temperature 


28 

28 

55 

51 

25 

500 mb height 


15 

18 

22. 

22 

20 
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over all regions analysed, 

!Ehe magniinde oJl the impact of the SSI anomalies on 
the model atmosphere over the ‘northern Hemisphere is in- 
dicated oy the rms differences and S1 comparison scores 
he’fcween the A and M simulations, shorn in Table 5 for 
each month. Compared with the simulation errors in 
Table 4|t the impact scores in Table 3 are virtually neg- 
ligible. The small S1 comparison scores for the 500 mb 
height field suggest that the impact of the SST anoma- 
lies on the model simulations is indeed trivial at that 
level. However, the influence of anomalous ocean tem- 
peratures does appear to be greater at lov^er levels of 
the model, and all the impact scores in Table 5 are above 
the “noise level” of the model, which is discussed below. 
Over certain smaller regions, the effect of SST anomsilies 
on the sea-level pressure pattern is, in fact, quite 
large. I'or example, in January 1977 the S1 comparison 
score over the United States between simulations A and M 
is 75 for sea-level pressure, indicating a considerable 
alteration of that field. On the other hand, the impact 
S1 score for sea-level pressure over the United States 
in February 1977 is only 48; and in neither month v;as 
the simulation of the sea-level pressure field improved 
over the United States by the use of observed SSTs, 

In general,, it appears that, although the use of 
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observed SSU}s i:n. place of climatoloslcal vaXues bad ao 
besieficial effect on tbe sirntLlations, tlie model is not 
insensitive to SS!I? anomalies. IKie model ^s response is 
negligible at upper levels, but not at sea level, where 
appreciable Iterations (not necessarily realistic) of 
the pressiore pattern may be generated both locall3!‘ and 
remoteli^ by SS2! anomalies, Ihe character of the local 
response is a negative correlation between SS2 anomaly 
and sea-level pressure change at the same point; i.e., 
cold (v;arm) anomalies generate higher (lov:er) sea-level 
pressures. Ihis relation is neither miiversal nor 
linear, but it is found with few exceptions. Por ex- 
ample, in a sample of .15 warm anomaly and 21 cold anomaly 
gridpoints with absolute magnitudes greater than 2° C in 
January 1977# only four of the 36 points did not e:chibit 
this behavior. Ihe m,ean values of the SSI anomaly and 
the Qo-located sea-level pressure difference between 
simulations A and M for the sample were + 5*0° 0 and 
- 2.3 mb for the warm anomaly points and - 3*2® C and 
+ 2.2 mb for the cold points. A synoptic effect of the 
SS5P anomaly field found in all months studied is a 
general lov^ering of pressures in the ITorth Atlantic and 
an elevation of pressures in the eastern North Pacific, 
but v/i’lh no improvement in the quality of the simulations, 
fhe direct, local, thermal effect of the S3T anomalies 
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be seen in tbe S 50 mb temperaijiare diiierences over 
the oceans between the A and H simulations* Por the 
sample above oi 15 warm and 21 cold anomaly gridpoints 
in d'anuary 1977» all bnt five points e:x)iibit a positive 
cornelation between the SS9} anomalies and the S50 mb 
temperature differences, with meaii values + 2,5° C and 
1*1*^ 0 at the 850 mb level corresponding to mean SSl^ 
anomalies of + 5-0*^ Q and 3*2° G, respectively 3}he 
large positive SSf anomalies, which generate strong up- 
ward heat fluxes from the sea surface, are generally 
associated with large positive temperature differences 
between A and I'l at the S5Q mb level, vmereas the stabi- 
lizing large negative SSf anomalies do not produce cor- 
respondingly large cooling effects on the model atmos- 
phere, n?he local decrease (increase) of sea-level pres- 
sure over warm (cold) SSD? anomalies is a hydx'ostatic 
conseq[uence of the local thermal effect. 

Idle total reaction of the model atmosphere to SSf 
anomalies is, of course, a complex non-linear combination 
of local and remote effects, with the latter ofteai domi- 
nating# (A more detailed synoptic analysis of the re- 
sponse of the model atmosphere to anomalous ocean tempera- 
tures will be presented in a separate publication.) 
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s - ■ 

Res-ponse to OIxaiiRes in InitiaXiga~feion. 

Model simulations of tiie atmospliere start with, ar- 
bitrarily selected initial conditions, IPor example, it 
is convenient to initialize a monthly' simulation run with 
data for. the first day of a calendar month. Such initial 
data may be viewed as a random sample tahen from the 
fluctuating state of a system whose characteristic time 
scale is 'of the order of days* How sensitive is a long- 
term mean simulation to this arbitrary choice of initial- 
ization? Hov; different would monthly mean simulations be 
if started with data one or two days earlier, or later? 
large differences, indicating an excessive sensitivity to 
the initial state, could ca^t serious doubt on the credi- 

hility of the model output. Jbi experiment was, therefore, 

% 

carried out with the climate model to determine its sensi- 

s 

tivity to systematic changes in large-scale initial con- 
ditions, such as might result from shifts of one or tv/o 
days in initialization time. 

In this experiment, the simulation runs for October 
1976 and January 1977 were repeated with the initializa- 
tion time shifted to 00 G-HT on days 2 and 5, respectively. 
[Che mean simulation from each run was then compared v/ith 
the mean for the corresponding period to the end of the 

g 

Adapted from a master's thesis submitted to Ihe City 
College by Robert Klugraan... 
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mojitlb. (mqn.1?l\^iTQiitx'Us-ono doy ?md mcmtli-mimis-wo days, xe-* 
spaatiYGl*Y) gauarated "by idia 'baaiQ simulation nun (tlia 
«cbntrQl«5 inom dajr i data^ QJlis principal xesi^lts oi 
tliis expanimant ana presented in Salle 6, in \'jnicn are 
snom tJie rms difierences (l)l and DS) 'between eacdi 
«sliiited iiidtialiaation*’ run and tie ''control*’, as well 
as tJio rms errors (Kl ojid 33S) of oaol run, aiid in Sable 
7, stowing the corresponding S1 comparison and still 
scores f or tte same simulations, for ’te Horttern Hemi-' 
sptej^e only*: 

Comparing 121 ajid 'B2. in Sables 6 and 7 with eaolx 
otter and witli M (tiie ’’control”) in lUables S and 5, it 
can be seen, first of all, that shifting initialisation 
time tad no significojat effeot on eitter 'bi^e rms errors 
or SI still scores of tiie montlily mean simulations, ill-' 
though tte verification periods of tte three simulations 
differ by one or two days, 'their rms errors and S1 s3:ill 
scores over the Hoi'thex'n Hemispliere are almost tl^e seiiie 
for corresponding mon-biis and variables, (Complete x^e- 
suits, not slxo^m, indicate that the sajjie is true fo,r all 
tlxe regions listed in fables 2 and 5*) 

Of ooxtrsQ, tile similarity of erx'ox^ scores does not 
necessarily mean tliab shifting initialisation time had no 
effect on tlie simulated synoptio fields, fhe rms diffex’- 
ences and Si comparison scores (hi and P2) in fables 6 
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DJable 6. Root-mean-sauare (rms) differences between 

I i' ' " '■■ ■* * 

ffiontbly mean simi’uLations initialised on day 2 
versus day 1 (Bl) and day 3 versus day 1 (B2) 
compared with corresponding rms errors for 
day 2 (El) and day 5 (jS2) initialisations for 
October 1976 and Janxxary 1977 over tbe Northern 
Kemispbere. (S^ee tex:t for details.) 

Sea-level 850 mb _ 500 mb 

Rressura (mb) gemperat\ire C 0) Heigjit (m) 




Oct. 176 

dan. '77 

Oct. '76 

Jan. '77 

0ct.»76 

Jan. '■ 

El 

(difference) 

1.4 

2.2 

1.5 

1.6 

19 

18 

El 

(err OX') 

5.1 

7.1 

4.9 

5.8 

99 

154 

1)2 

(difference) 

1.9 

5*1 

1.8 

2.2 

25 

53 

E2 

(error) 

5.1 

7.2 

4.S 

5.5 

95 

151 



gable 7« S1 comparison scores for moniilily moan simu- 
lations initialiaed on day 2 yorsus day 1 (1)1 ) 
and day 5 versus day 1 (D2), comparod with cor- 
responding S1 shill scores for day 2 (El) and 
day (E2) initialisations for Octobor 1976 
and January 1977 ovei^ the ITorthorn Uemisphore* 
(See text for details.) 


SQa*^level 

850 mb 

500 mb 

Pressure 

gomoorature 

HeiRht 

Oct. "76 Jan. "77 

Oct. *76 Jan. "77 

Oct. "76 Jan. "77 


1)1 (comparison) 

32 

40 

23 

27 

19 

19 

El' (shill) 

92 

87 

62 

65 

50 

67 

P2 (comparison) 

43 

51 

55 

33 

25 

27 

E2 (shill) 

91 

86 

59 

64 . 

50 

69 


and 7 ixidica1>e quantitaijively the magnitude of this effect. 
It is apparent that the differences due to shifting ini- 
tialisation time are smaller than the simulation errors* 
HoweTer, the magnitude of the effect is not entirely neg- 
ligible* Furthermore, it is greater for a two-day (D2) 
than for a one-day (Dl) shift, 33ie effect at the 500 mb 
level is almost trivial, especially in terms of the S1 

4 

comparison scores; but the same is not true of the other 
tv7o fields. !Ehe initialization shift effect is also 
smaller than the simulation error of climatology (C in 
fables 2 and 5). However, the rms difference betv/een 
850 mb temperature simulations with a two-day shift in 
initialization (l)2 in fable 6) is only slightly smaller 
than the corresponding simulation error of climatology 
(Q in fable 2), 

In general, shifts of up to two days in initializa- 
tion time produce small, but not trivial, changes in the 
monthly mean synoptic fields generated by the climate 
model, fhe magnitude of this effect represented by the 
values of hi and 1)2 in fables 6 and 7 indicates the in- 
herent minimal error of the monthly mean simulations as- 
sociated with the arbitrary choice of initialization time, 

fhe use of instantaneous initial conditions for the 
.generation of long-term mean ‘simulations raises other 
questions regarding the role of small-scale and 
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short-lived compoxieiits in the initial analysis* Do such 
transient features of the initial state adversely affect 
the mean, simulations, and would the use of a time filter 
applied to the initial data improve the result? 

An experiment was carried out with the October 1976 
and January 1977 data in which the five 12-hourly ob- 
served fields for the first 2-^ days of each month were 
first averaged arithmetically. 'Bie model was tlien ini- 
tialised at 00 &MI2 on day 2 with these time-avex'aged 
data and run to the end of the calendar month* Hie re- 
sulting monthly (minus one day) mean simulations were 
then compa3.'ed with the observed means for the same periods, 
with the results shown as M in Table 8 . j\1so shov'ni for 
comparison in the table are the corresponding error 
statistics, M, for the basic simulation run, fi-om Tables 
2 and 5. Although the results are shown only for the 
northern Hemisphere in Table 8, they are essentially the 
same for all the regions listed in Tables 2 and 5. 

It is obvious from Table 8 that the use of time- 
averaged initial conditions did not result in any con- 
sistent improvement in the monthly mean simulations. The 
rms errors and SI skill scores for K a3.'e essentially the 
same as for M, and in some cases slightly worse, (Al- 
though they are not shorn* the synoptic maps as well as 
the rms differences and 31 comparinon scores for h vs, M 


57 



5?able S. Root-mean- square (rius) errors and S1 skill 


scox-es of monthly meaia simulations initialised 
with time-averaged initial data, H, compared 
with those from day 1 initial conditions, M, 
(ii'om lahles 2 and 3) for Octoher 1976 and 
January 1977 over the Northern Hemisphere. 

(See tej:t for details.) 


X'lns error 

0ct.»76 Jan. ’77 

Sea-lovel 

Pressure 

Oct. ’76 Jal^.’77 

850 mb 
temperature 

Oct. ’76 
500 mb i 

Jon. ’77 
lei?rht 


mb 

°0 

m 


H 

4.9 

7.5 

4.9 

6,0 

102 

156 

M 

5.5 

7.2 

5.1 

5.9 

104 

131 

S1 score 







I 

91 

89 

62 

65 

55 

68 

H 

95 

89 

66 

67 

52 

67 


5a 



also indicate no significant differences between tlie two 
simulations.) !Eius, no benefit was gained from the nse 
of time-averaged initial data. 

njhe poor quality of monthly mean model simulations 
is a consequence of the rapid decay of predictability, 

An experiment was carried out to determine if this decay 
could be retarded by averaging the output of the first 
five days of the basic computation, then re-initialising 
the model at the middle of the five-day period with these 
time-averaged output data for the remainder of the calcu- 
lation. Montlily simulations for October 1976 and January 
1977 were computed by this method, -v/ith the first tliree 
days of the mean talcen- from the basic simulation i^un and 
the remainder from the r e-ini tialiaed output, (The veri- 
fications of these simulations, together with those of» 
the basic experiment (from Tables 2 and 3)» are show in 
Table 9 for the Northern Hemisphere. Also shovoi in the 
table are the rms differences and S1 comparison scores 
between the monthly mean simulations computed by the two 
methods, where M* represents the result of the re- 
initialisation procedure, and M, as before, the basic 
simulation.- 

It is apparent from Table 9 that the re-initialisation 
method produced no discernible change in the monthly mean 
simulations. Not only are the rms errors and S1 skill 
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gable 9 


^ Hoot-mean- square (rms) errors and S1 skill 
scores of monthly mean simulations for October 
1976 and January 1977 over the Northern Hemi- 
sphere by the re-initialization method, H*, • 
compared vmth the basic simulation, M. ilso 
sho-wn are rms differences and S1 comparison 
scores betvjeen M* and M. (See text for details.) 


Oct. *76 Jan. *77 

Sea-level 

Pressure 


rms mb 


M* error 

5.1 

7.4 

M error 

5.3 

7.2 

(M* - M) difference 

0.9 

1-5 

S1 



M* skill 

93 

90 

M skill 

95 

89 

(M* " M) comparison 

22 

26 


Pet. *76 Jan. *77 Oct. *76 Jan. *77 


850 mb 
temperature 

500 mb 

Keisiht 


°0' 

m 

5.1 

6.0 

107 

135 

5.1 

5.9 

104 

131 

0.9 

1.5 

11 

% 

19 

65 

66 

52 

65 

66 

67 

52 

67 

16 

18 

10 

14 
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scores almost identical for the two methods, hut the dif~ 
ferences between the two simulations, as indicated by the 
rms differences and SI comparison scores are practically 
negligible. (Similar results, not shown, were found for 
all seven regions.) ^Ehus, neither time-averaging of ini- 
tial condi'" vons nor re-initialization with time-averages 
of early model output appears to result in any improve- 
ment in the quality of the monthly mean simulations. 
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7 ' 

Hoise Level Estimation' 

Inherent imcertainties exist in all model simulations 
due to unavoida'ble random errors in the initial data, 

!Sie influence of such errors on long-term atmospheric 
simulations has been evaluated in "noise level” experi- 
ments (e.g,, Chervin and Schneider, 1976 a, b; Spar et 
al,, 1978) in which models were run repeatedly with the 
same synoptic-scale initial conditions, but with the ini- 
tial analysis contaminated for each simulation by a dif- 
ferent small-scale random perturbation field, !Ehe dis- 
persion of the resulting simulations represents the noise 
of the model, which should be weak if it is not to mask 
signals due to prescribed changes in the model, 

A noise level calculation was carried out with the 
coars e-mesh climate model to determine the inherent un- 
certainty of the monthly mean simulations. In this ex- 
periment, the initial analysis for 00 GMI 1 October 1976 
was perturbed four times by means of a random number 
generator which geographically distributes Gaussian 
error sets with specified global rms magnitudes, Ihe 
global rms errors added to the initial data in this ex- 
periment were 5 Jes for the zonal and meridional wind 
components and G for the temperatures at all nine 

'^Adapted from a master *s thesis submitted to ahe City 
College by Jesus J, Kotario, 
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sigma leyels of the models and 5 mb for the surface pres- 
sures, I’our prediction, runs fr'om the perturbed initial 
states were then executed, resulting in a total set of 
five monthly mean simulations , ' including the basic simu- 
lation which serves as a control, 

ITom a visual comparison of the five fields of sea- 
level pressure, 85 O mb temperature, and 500 mb height 
(not reproduced), there appear to be no discernible 
synoptic differences among the simulations, The rms er- 
rors and S1 skill scores (not shown) of the four perturba- 
tion simulations differ only slightly from those of the 
control shown in Tables 2 and 5« The similarity of the 
five fields is indicated quantitatively by the small rms 
differences and S1 comparison scores between pairs of 
simulations over the northern Hemisphere in Table 10, ' 
where the results of the first perturbation, ]?1, are com- 
pared with each of the other four runs. Here it can be 
seen that the averages of the rms differences over the 
hemisphere are approximately 1 mb, 1° 0, and 10 m, re- 
spectively, for sea-level pressiire, 850 mb temperature, 
and 500 mb height. These values are less than half as 
large as those found by Spar et al, (197S) in the cor- 
responding experiment v/ith the 4x5 degree GISS model, 
probably reflecting the influence of the. coarser resolu- 
tion. Tlie averages of the S1 comparison scores in 
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IPahle 


Sea-level 

Eress'irpe 

8^0 mb 
Temperaijure 

500 mb 
Height 


10. Hoot-mearx-rspuare (arms) differences and S1 

comparison scores over the Northern Hemisphere 
between monthly mean simulations for October 
197^- M denotes control simulation. PI, P2, 
P5, and P4 represent simulations from four 
different random perturbations of the control 
initial conditions. 


PI - M 


rms difference (mb) 0,8 

S1 comparison 25 

rms difference ('^O) 0.9 

S1 comparison 16 

rms difference (m) 9.1 

S1 comparison 10 


PI - P2 PI - P5 PI - P4 

1,0 1.1 1,5 

26 27 29 

0.9 0.9 1*0 

18 17 19 

% 

9*0 9.4 11,8 

10 11 14 
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lU^ble 10 are approximately 26, IS, and 11 for the tliree 
fields, again indicating virtually no differences among 
the hemispheric patterns due to the initial perturbations. 
It is apparent from lable 10 that, in genera" the in<" 
fluence of random initial state errors on the monthly 
mean simulations is negligible compared v/ith the simula- 
tion errors. Similar results were found for all regions 
analysed* 

idle geographical distribution of model noise level 
may be repi*esented on global maps of standard deviation 
of the simulated variables (Chervin and Schneider, 1976 b). 
Ihese were computed at each gridpoint in the present ex- 
periment from the five, simulations (with four degrees of 
freedom) for sea-level pressure, 850 mb temperature, and 
500 mb height, and are slio\m in I'igures 5, 6 and 7, re- 
spectively, in digital form. Longitudes are labelled at 
top and bottom, negative for west and positive for east, 
with the Western Hemisphere on the left. Latitudes are 
sliomi in the second column from the right, positive for 
north and negative for south* JSonal mean values of the 

standard deviations are tabulated in the first column on 

—1 —1 

the right, and the units on the maps are 10 mb, 10 
degi'ees 0, and m. (Also shoMi on the maps are the 24 x 56 
gridpoint nuiabers . ) 

5tom the Bonal mean standard deviations tabulated 
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5, Standard deviations of five simulations of the global sea-level pressure field for 
October 1,976. 

Units': 10 ^ mb. 
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Figure 6. Standard deviations of five simulations of the global 850 mb temperature field f or . 
October 1976. 
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Units: 10 degrees C. 
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Figure 7. Standard deviations of five simulations of the global 500 mb geopotential height field 
for October 1976. 

Units : meters . 


in I'igures 5, 6 and it is apparent tlaat, in general, 
tlie noise level is negligibly small near the Equator 
(with zonal mean standard deviations of 0.2 mb, 0,4^ C, 
and 1 m) and greatest in high latitudes, Ihe largest 
point standard deviations of sea-level pressure are found 
in the Antarotio (3*6 mb) and near Iceland (5*0 mb), with 
maximum zonal means of 1.8 mb near the South Pole and 
1,2 mb at 74° Por the 850 mb temperature simulations, • 
the largest point standard deviations are found again in 
the Antarctic (6,2° C) and also over Greenland (5.5° C), 
with maximum zonal means of 2.0° 0 at 82° S and 1.0° C at 
74° N, At the 500 mb level, maximum point standard de- 
viations appear in the Antarctic (20 m). South Atlantic 
(19 m), Eorth Atlantic (21 m), and over eastern Siberia 
(24 m), with maximum zonal means of about 10 m north of 
51° N and south of 42° S, 

!Che increase in noise level from low values in the 
tropics to higher values in high latitudes is consistent 
with the results of Chervin and Schneider (1976 b) v;ith 
the ECAP* model and Spar et al. (1978) v;ith the GISS model, 
However, there is considerable zonal variation of the 
standard deviations as v;ell in Figures 5, 6 and 7, even 
in the sub-tropics, v/here larger values generally appear 
over the continents (Australia, Africa, South America) 
than over the oceans. lEhis is particularly true of the 
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850 jti’b ’tjemperatures, indicaiiins the apparently greater 
eeneltiYity of the model temperatwoe aaXoulations to 
.initial oonditions oyer land than over v/ater. The noise 
level also appears greater over the Antarctic and G-reen- 
land ice* possibly reflecting modeling defects associated 
with topograpliy as well as the physics of the ice surface- 
In general, it may be concluded that, v/hile the in-- 
fluence of random initial state errors on the monthly 
mean simulations appears to he insignificant when com- 
pared with the large-scale simulation errors, the geo- 
graphical distribution of noise level, particularly its 
variation with latitude, should be considered, as suggested 
by Ohervin and Schneider (1976 b), in the evaluation of 
any climate change experiments with the model, Hov/ever, 
the relatively trivial influence of random initial state 
errors compared with the simulation errors indicates both 
stable behavior of the model and the need for further 
model improvement. 



Ooiioluaious 


Mosatlily meaii siniulaijions computed fx'om initial data 
at the beginning of each of the five montlis, October 1976 
through ]?ebruary 1977 » with a coarsG“resolution climate 
model developed by Hai^son mid his colleagues at OISS have 
revealed certain limitations of the model as a long«range 
prediction system, despite its realistic simulation of 
largG'-scale meridional structure. However, as the model 
is still in the process of development and modification, 
any conclusions regarding its performance should be con- 
sidered only tentative. 

3}he model simulations exhibit a cold bias in the 850 
mb temperatures at all latitudes, with a corresponding 
negative bias in the geopotential heigi^ts of the 500 mb 
surface, relative to climatology as well as to the ob- 
served monthly mean state. For most of the regions, and 
months analysed, the inns errors and SI skill scoi'gs of 
the monthly simulations of sea-level pressure, 850 mb 
temperature, ai^d 500 mb height are larger than the eoi’- 
responding values for climatology, indicating that, at 
this stage, the model does not realistically simulate 
anomalous monthly mean synoptic patterns from given ini- 
tial conditions, 

Mion observ'^d monthly nv(;rage soa-surface tempera- 
tiircs (SSTs) are substitutea for the climatological 



SS3!s used in the .first simulations, the results are not 
improved. At upper levels, represented by the 500 mb 
surface, the impact of SSff anomalies on the model simu- 
lations is found to be negligible, lower levels of the 
model atmosphere, on the other hand, are more sensitive 
to alterations in ocean temperatures. Thus, for example, 
higher (lower) sea-level pressures are generated locally 
by tlae model over cold (warm) SST anomalies, and more 

4 

complex remote effects on the sea-level pressure field 
are fomid also. However, the impact of the SST data on 
the monthly mean simulations of sea-level pressux'e, 850 
mb temperature, and 500 mb height fields is, in general, 
not beneficial, and the results of the experiment lend 
no support to the hypothesis that SST anomalies are re- 
sponsible for atmospheric anomalies. Admittedly, this' 
resxilt may be more indicative of model deficiencies than 
of the behavior of the real atmosphere. 

Small changes in the monthly mean synoptic patterns 
generated by the model are found to result from shifts 
of one and two days in the initial data, indicating an 
inlierent minimal error level associated with the arbi- 
trary choice of initialisation time. The monthly mean 
simulations are not improved by the use of either time- 
avez'aged initial data, smoothed over a 2i'’-day period, or 
re-iiiitialiaation, before completion of the monthly run, 
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with time-averaged model output smoothed over the first 
five forecast days. 

SChe noise level of the coars e-mesh model, as deter- 
mined from a multiple random perturbation experiment, 
appears to be even lover than that of the 4x5 degree 
GISS model. At this stage, the general effect of random 
initial state errors on the monthly mean simulations is 
negligible compared with the large-scale simulation er- 
rors, indicating stable behavior of the model but also 
a need for model improvement. Average rms differencet 
betveen perturbed simulations over the Northern Hemisphere 
are only about 1 mb, 1° C, a:id 1 m for sea-level pressure, 
850 mb temperature, aiid 500 rab height, respectively. 
However, the global distribution of the noise level, 
represented by maps of the standard deviations of the 
five simulations genex'ated in the perturbation experiment, 
reveals marlced spatial differences, with minimum values 
in the tropics and maxima in high latitudes, as well as 
zonal variations which suggest greater sensitivity to 
initial errors over land and ice than over open ocean. 
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yiKures 


1, Monthly mean sea-level pressiires over the northwest 
quadrant. !Eop; observed; bottom; simulated. Units: 
mb “ 1000, 4 mb isobars, (a) October 1976 through 

(e) 3?ebruary 1977* 

2* Monthly mean 850 mb temperatures over the northv/est 
quadrant. 0?op: observed; bottom: simulated. Units: 
degrees 0. 5^ G isotherms, (a) October 1976 tlarough 

(e) ]?ebruary 1977. 

5. Monthly mean 500 mb geopotential heights over the 
nortlwest quadrant, Top; observed; bottom; simu- 
lated. Units; decameters - 500, 100 m contours, 

(a) October 1976 tlirough (e) February 1977* 

4. Monthly mean sea-suriace temperature anomalies 

(observed - climatology) over the globe in degrees 
0. (a) October 1976 through (e) February 1977* 

5. Standard deviations of five simulations of the global 
sea-level pressure field for October 1976, 

Units; 10"*^ mb, 

6. Standard deviations of five simulations of the global 
850 mb temperature field for October 1976, 

Units; 10*’’^ degrees 0, 
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standard deviations of five simulations of the global 
500 mb geopotential height -field for October 1976. 
Units : m. 



